Modified Coulomb law in a strongly magnetized vacuum 
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We study electric potential of a charge placed in a strong magnetic field B 3> Bo — 4.4- 10 13 G, as 
modified by the vacuum polarization. In such field the electron Larmour radius is much less than its 
Compton length. At the Larmour distances a scaling law occurs, with the potential determined by a 
magnetic-field-independent function. The scaling regime implies short-range interaction, expressed 
by Yukawa law. The electromagnetic interaction regains its long-range character at distances larger 
than the Compton length, the potential decreasing across B faster than along. Correction to the 
nonrelativistic ground-state energy of a hydrogenlike atom is found. In the limit B = oo, the 
modified potential becomes the Dirac <5-function plus a regular background. With this potential the 
ground-state energy is finite - the best pronounced effect of the vacuum polarization. 

PACS numbers: 03.50, 12.20.-m, 77.22.Ch, 97.60.Gb 



There is now compelling evidence that many compact 
astronomical objects (soft gamma-ray repeaters, anoma- 
lous X-ray pulsars, and some radio pulsars) identified 
with neutron stars have surface magnetic fields as high 
as ~ 10 14 — 10 15 G More strong magnetic fields 
(B ~ 10 16 — 10 17 G) are predicted to exist at the surface 
of cosmological gamma-ray bursters if they are rotation- 
powered neutron stars similar to radio pulsars 0]. All 
these fields, however, are much smaller than the maxi- 
mum value inherent in quantum electrodynamics Q . 

Vacuum in an external magnetic field B behaves as 
an anisotropic dielectric medium with spatial and fre- 
quency dispersion, (e.g., [1]). These properties may be- 
come important, provided that the field strength achieves 
the characteristic value Bq = m 2 /e ~ 4.4 x 10 13 G, where 
m is the electron mass and e is its charge. [Henceforth, 
we set h = c = 1 and refer to the Heaviside-Lorentz sys- 
tem of units.] Although much work has been devoted to 
study of electromagnetic wave propagation in the magne- 
tized vacuum, problems of electro- and magneto-statics 
in this medium did not attract sufficient attention, save 
Refs. [f|, where corrections to the Coulomb law were 
found when these arc small: for B/B -C 1 in d, or at 
large distances from the source for 1 <C B/Bq <C 37ra -1 
in @ (a = e 2 /47r = 1/137). 

In this Letter, we find that for sufficiently large b = 
B/Bq 3> 1 the electric field produced by a pointlikc 
charge at rest may be significantly modified by the vac- 
uum polarization, the modification being determined by 
the characteristic factor ab. The modified Coulomb po- 
tential in the close vicinity of its charge, characterized 
by the Larmour length L-q = (e£?)~ 1//2 = (1/myb), goes 
steeper than the standard one, following a Yukawa law, 
whereas it obeys a long-range " anisotropic Coulomb law" 
far from the source, at distances characterized by the 



electron Compton length m 



^> Lb- (Details of 



the corresponding derivations can be found in the accom- 
panying preprint 0]). The short-range part of the mod- 
ified potential tends to the Dirac (5-function in the limit 



b — ► oo. The modification of the Coulomb law should af- 
fect, first of all, the field of an atomic nucleus, placed 
in a magnetic field. We determine the corresponding 
correction to the lowest energy level of a hydrogenlike 
atom. Unlike the famous result of Ref. Q , referred to in 
many speculations on behavior of matter on the surface 
of strongly magnetized neutron stars (e.g., the review 
Q and references therein), we find that the (vacuum- 
polarization-corrected) ground state energy remains fi- 
nite in the limit of infinite magnetic field. 

Let the constant and homogeneous magnetic field B 
be directed along axis 3 in the frame where the pointlike 
charge q is at rest in the origin x = {x\, x%, x^\=Q, and 
no external electric field exists. By using the tensor de- 
composition of the photon propagator over eigenmodes in 
a magnetic field 13, HI it is straightforward to show Q 
that electrostatic potential A$ produced by this charge 
has the form 

exp(— ikx) d 3 fc 



(27T) 



(i) 



= 0. The 



while its vector potential is zero, j4i ) 2,3(x) 
static charge gives rise to electric field only, as it might 
be expected. Here k 2 is one (out of three) eigenvalue(s) 



(a) 



of the polarization operator H^: 

1,2,3, fi, v = 0,1,2,3. The eigenvectors are 4- 
potentials of the eigenmodes. The eigenvalues depend 
on two combinations of the photon momentum compo- 
nents k% — and k\ = k\ + k\ taken at zero frequency 
ko = 0. Eq. ([1]) is approximation-independent and axial 
symmetric: A (x) — A (\x3\, x±), x± — \J x\ + x\. 

Eq. ((T|) indicates that only mode-2 photons mediate 
electrostatic interaction. This fact may be better under- 
stood, if we examine electric and magnetic fields intrinsic 
to the virtual (off-shell) photons of this mode, obtained 
from its 4-vector potential \> v = (kz ko) v . It may be 
seen from explicit representation for these fields ll| that 
in the limit ko — the magnetic field of mode-2 photon 
disappears, while its electric field is collinear with k, i.e., 



2 



it becomes a purely longitudinal virtual photon. Virtual 
photons of other modes, in the static case, are carriers 
of stationary magnetic fields. For instance, mode-1 pho- 
tons are responsible for the field produced by a constant 
current, flowing parallel to the external magnetic field. 

In the asymptotic regime of high magnetic field, eB 3> 
fc| and B ^> m 2 /e = B , with the accuracy to terms 
that only grow with B as its logarithm or slower, the 
eigenvalue K9, calculated within the one- loop approxi- 
mation in [l(J, U I , acquires the form 0, Hit ] 



2abm 



■ exp 



2m 2 b J V Am 2 J ' 



T{y) = V 



(1 - v 2 )dy 

1 + 2/(1 — 77 2 ) ' 



(2) 



(3) 



Note the properties: T(y -> 0) ~ 2y/3, T(oo) = 1. 

The deviation of the potential iJTJ from the standard 
Coulomb potential Aq (x) = q/{4jry/Tj_ + x 3 ) 



is 



AA (x) = ^(x)-A (x) 

/ exp(-ik 3 x 3 ) 

> V 



Jo{k±x±)dk 2 



exp(-ik 3 x 3 ) 



hi 



fe 3 



k ~\~ k^ K2{k^-} kj^) 



dk 3 . (4) 



Here Jo is the Bessel function of order zero. This integral 
defines AA (x) as a finite function of the coordinates in 
the origin, unlike A and Aq. As the integration variable 
k 3 in (fj| approaches the large values ±y/eB, one has 
k-2 <C k 3 . Hence, no essential contribution comes from 



the integration over the region \ks\ > y/eB, wherein eq. 
d2|) is not valid. The expression to be substituted for it 
there is even smaller (since it does not contain the large 
factor b) [l3| • The leading terms of the expansion of ([TJ) 
near the origin x 3 — x±_ = are 

A (x) ~ ± ( ±r - 2mc], C = — AA (0)>0, (5) 
47T \|x| / qm 

where C is a constant depending on the magnetic field. 

More exactly the singular behavior near the origin will 
be presented if we note that it is provided by the inte- 
gration over large k 3 (and k±) in JT]), where we may set 
T{k 3 /Am 2 ) = T{oo) = 1 and perform ^-integration in 
Eq.([T]) by calculating residues. Then 



^o(x) 



Aq(x) 
Lb 



J {k±x±)kj 



exp 



|s 3 |W*i + (2a/7r)exp [-k 2 j2 



k 2 L + {2a/Tr) exp(-fc^/2) 



dfc 



(6) 



Here Ao(x) is a dimensionless, external-field- independent 
function of the arguments x 3 = x 3 /Lb, %± = x±/Lb, 



and the integration variables k 3 = k 3 L^ 1 k± = /cj_Lb are 
used. We refer to Eq. (J6J) as establishing a scaling regime 
that describes the potential measured in inverse Larmour 
units as a universal function of coordinates measured in 
Larmour units. It holds for |x| <C (2m) _1 . 

The simple representation ([6]) can be further simplified 
if x 3 or x± are large in the Larmour scale: \x 3 \ 1, or 
^> I (but remain small in the Compton scale). In 
this case the integration in ^ is restricted to the do- 
main fc^ < 1 where the exponential exp(— k 2 _/2) should 
be taken as unity. Then ^ is reduced to the isotropic 
Yukawa law 



A U) 



exp {~(2a/TT) 1 / 2 ^/x 2 L +xfj 



4itLt 



(7) 



This can be established by tracing © back to ([T]) with 



k 2 (oo, 0) 



2a 



2ab 



M 



(8) 



substituted for —K2(k 3 ,k 2 J in the denominator. Here 
M is the "effective photon mass" noted in Ref. [IBj]. 
The Yukawa law ([7]) establishes the short-range charac- 
ter of the static electromagnetic forces in the Larmour 
scale. Stress, however, that the genuine photon mass 
understood as its rest energy is always strictly equal to 
zero as a consequence of the gauge invariance reflected 
in the approximation- independent relation K a (0, 0) = 
respected by @. Hence, the potential, produced by a 
static charge, should be long-range for sufficiently large 
distances. This is the case, indeed. One can see by in- 
specting the curves of Fig.l, computed using Eq. (fT|) 
for x± — 0, that the scaling regime ©, or (J7J docs 
fail for sufficiently large values of x 3 : ~ 0.1/ (2m) — 
50L B for b = 10 6 , - 0.2/(2to) = 33L B for b = 10 5 , 
~ 0.3/(2to) = 15Lb for b = 10 4 - the larger, the smaller 
the field [for these distances Eqs. © and (J7J are al- 
ready the same]. Starting with these values, the poten- 
tial curves approach their envelope, that can be fitted as 
A (x 3 ,0) ~ ( 9 /47r)1.41/(x 3 + 1.04/2m), unlike the scal- 
ing curves ([6]), and ([Jj that tend fast to zero. It is in this 
place that the abruptly falling - short-range - potential 
turns into a slowly decreasing - along the envelope curve 
- long-range potential. An analogous change from the 
short- to long-range behavior is observed in Fig. 2, where 
the electron energy in the field |T]) is plotted against the 
transverse distance from the charge at x 3 =0. At small 
distances all the dashed curves in Figs. Q] and [2] approach 
the thick solid Coulomb curve in accord with ([5|) . 

The scaling regime (]6]), (0, depends on the fact that 
the eigenvalue ([2]) grows linearly with the magnetic field. 
If this linearity, supposedly, retains in higher-loop ap- 
proximations, one may conjecture that the calculations 
of the latter may be reduced to finding a"-corrections 
to the mass ([5]). Potential ^ implies the suppression of 



3 



electrostatic force by the linearly growing term in the de- 
nominator of ([1]) at large distances in the Larmour scale 
from the charge, but not close to the charge, where it has 
the same singularity as the Coulomb law. 

For larger distances x 3 , x±, not small in the Compton 
scale, one may expect that only integration over small 
k 3 , k± is important in ([T]) (for a more thorough analysis 
see 0)- In this limit © behaves as k 2 (A:|,0) ~ -§ffc 3 - 
With this substitution the integral in fl} is 



A (x 3 ,x±) 



(J 



47r y/(x' ± ) 2 + x 



(9) 



1 /2 

where x' ± = (3x±, f3 = (1 + a6/37r) , x ± > x±. For 
small ab/i-K, Eq. © coincides with the result of [fj]. 

Eq. © is an "anisotropic Coulomb law", accord- 
ing to which the attraction force decreases with dis- 
tance from the source along the transverse direction 
faster than along the magnetic field, but remains long- 
range. In accord with ((9|), the curves Aq(x3,0) in 
Fig[T] all approach at large distances \x 3 \ the Coulomb 
law q/(47r|a; 3 |), whereas each curve A (0,x±) in Fig. 
[2] reaches at large x± the asymptote Aq(Q,x±) = 
q/(4:Trx' ± ) = Aq (0, x±)/[3, different for each field, i.e., the 
potential is anisotropic. Again, the same as for short dis- 
tances considered above, we face - now anisotropic - sup- 
pression of the Coulomb force due to the linearly grow- 
ing term in ([2]). The equipotential surface is an ellipsoid 
stretched along the magnetic field. The electric field of 
the charge E = — VAq(x 3 , x±) is a vector with the com- 
ponents (q/2n)(x 3 + /3 2 x 2 L )~ 3 / 2 (x 3 , f3 2 x±). It is not di- 
rected towards the charge, but makes an angle <f> with the 
radius-vector r, cos</> = 2 x 2 L )(x'l + (3 A x\)~ 1 / 2 (x 2 + 
x 2 ^)^ 1 / 2 . In the limit of infinite magnetic field, (3 = 00, 
the electric field of the point charge is directed normally 
to the axis 23. 

The regime ([9]) corresponds to the approximation, 
where only quadratic terms in powers of the photon mo- 
mentum are kept in «2. Within this scope the dielectric 
permeability of the vacuum is independent of the fre- 
quency, and the refractive index depends only upon the 
angle in the space [la |. 

Consider an impact of the Coulomb potential modifi- 
cation on the nonrelativistic hydrogenlike atom ground 
state energy. The wave function 'J (23) is subject to one- 
dimensional Schrodinger equation [8[ with respect to 23, 
valid in the region \x 3 \ > Lb- The transverse coordinate 
x±_ in the argument of the potential is replaced by the 
momentum of the transverse center-of-mass motion [13] 
and should be set equal to zero for the ground state. 

For the fields in the range 1 <C b <C 2-7r/a ~ 10 3 the 
correction (j4]) to the Coulomb potential may be treated 
as perturbation. Keeping only the first-order term in «2 
in Q we may derive the magnetic analog of the Uehling 
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FIG. 1: Electron potential energy —eAo in the field of a point 
charge q — Ze plotted in units 2aZm — Zx7.46 keV as 
a function of longitudinal coordinate X3 in Compton half- 
lengths (2m)" 1 at ii =0 for four values of magnetic field 
(from left to right): b = 10 6 , 10 5 , 10 4 , 10 3 , b = B/B , 
B = 4.4 x 10 13 G (dashed lines). Solid line is the fit 
— 1.4/(2mx3 + 1.04). Bold line is the Coulomb law — l/2mxa. 
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FIG. 2: The same as in Fig[T] in function of transverse co- 
ordinate x± at X3 — 0. The asymptotes at large x± values, 
— l/2mx'±, are presented by solid lines. 



potential [3]. At x±_ — it is 
AA (x 3 ,0) ~ ^ f' 2 e -( 2 ™l-3l/cos0) CQS 2 d0 (10) 



3vr 2 Jo 

This correction is of the order of ab, i.e., is much larger 
than a. Eq. (jlOjl implies C ~ abjVo. Calculating its ma- 
trix element with Loudon's Q wave functions we obtain 
the positive correction (q = Ze) 

2Z 2 a 3 bm, b / Inb \ . , 

In —j = Z 2 b — — Tlx 0.356 eV (11) 

37T 4a 2 V 8.454 / v ' 



to the negative Loudon-Elliott ground state energy 

Vb 

2a 



E 



-2Z 2 a 2 m\n 



(12) 



Eq. (|11[) may be used for Z < 6. For instance, for Z = 1 
and b = 200 it makes about 100 eV, while Eq ~ —2.5 
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keV. Although the proton has a finite size R ~ 10~ 13 cm, 
the Coulomb 1 / |a?3 1 part of the potential ([5]) remains the 
same within the definition range of the Schrodinger equa- 
tion, since R < Lb within the range of b considered in this 
paragraph (moreover, R <§; Lb). On the contrary, the 
vacuum-polarization part does depend on R, like it does 
in no- magnetic-field case [l8j]. However, the finite-size 
correction to (TTTj) is ~ (R/Lb) 2 to be neglected within 
the present scope of accuracy. 

The singularity 1 /|aj3 1 of the Coulomb potential in 
the origin is known to lead to the energy spectrum un- 
bounded from below 0: as this singularity is cut off at 
the Larmour length, Lb = (mVb)^ 1 , the ground state 
energy (fT2|) tends to — oo with the growth of the mag- 
netic field, when Lb — ► 0. This feature is cured by the 
vacuum polarization. As b — ► oo the region in Fig.l 
between the potential curve and the abscissa below the 
point —eV ~ — 1.4 x 2aZm, where it is crossed by the 
envelope, becomes infinitely deep and thin. The area 
S = (47r/(7) Jig A (x 3 , 0)dx3 calculated with expression 
([5]), where X3 is found from the equation A (x 3 ,0) = V, 
has a finite limit if and only if C grows proportionally 
to 1/Lb. This is the case: when calculated following Eq. 
@, C ~ 0.9594^/a6/27T. Thus, the dominating behavior 
of the modified Coulomb potential in the origin becomes 
the ^-function. Combining it with the fit for the envelope 
we come to the limiting form of the electron potential en- 
ergy at b — 00 (here |x3| does not exceed a few (2m) _1 ) 



eA (x 3 ,0) = -2aZ S <5(x 3 ) + 



1.4m 



2m\x 3 \ + 1.04 



(13) 



where S = ln( A /7r/2a/0.96) - 1 + 0.96 v /2o7^f = 1.79. 
A more rigorous calculation done with the use of eqs. 
(|6|) or for the short-range part of the potential in- 
stead of © results in S = 2.18 » -Ei(-(2a/ 7 r) 1 / 2 ), Ei 
is the exponential integral. With the l/|x|-singularity 
replaced by the <5-function, the ground energy level is 
certainly finite. Applying the formula for the ground 
energy — 2m(J^ A (x 3 , 0)dx3) 2 valid in a shallow well 
potential [h|, with eq. (Tl3|) integrated up to the point 
X3 = 2.6/2™, where the fitted curve crosses the Coulomb 
law, and the latter taken as the integrand for larger X3 up 
to the Bohr radius <zb = (ma)^ 1 , we estimate the finite 
limiting value for the ground energy as 



L, im = -2mZ 2 a 2 73. 



= -Z 2 x 4 keV. 



(14) 



The Loudon-Elliott energy ()12|) would overrun the limit- 
ing energy (fT4|) already for the magnetic field as large as 
b = 6600, when yet R -C Lb- The ground level reaches 
92% of its limiting value for b — 5 x 10 4 . After the mag- 
netic field reaches the value b = 1.5 x 10 5 , when R and Lb 
equalize, the Coulomb potential is cut off at the proton 
size, X3 = R. Setting Lb = R in (fT2|) we would get the 
minimum value for the Loudon-Elliott energy (Z = 1) to 
be —5.6 keV, which is essentially lower than (Tj"4")) . 



The significant modification of the Coulomb potential 
of an electric charge by the vacuum polarization in exter- 
nal constant magnetic field B ^> L?o, shown to eliminate 
the unboundedness from below of the nonrelativistic hy- 
drogen spectrum, is apt of having more implementation 
as far as other situations where electrostatic fields are 
important are concerned, e.g. properties of matter on 
surfaces of extremely magnetized neutron stars. 
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